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Spin-chain models have been widely studied in terms of quantum information processes, for instance 
for the faithful transmission of quantum states. Here, we investigate the limitations of mapping this 
process to an equivalent one through a bosonic chain. In particular, we keep in mind experimental 
implementations, which the progress in integrated waveguide circuits could make possible in the 
very near future. We consider the feasibility of exploiting the higher dimensionality of the Hilbert 
space of the chain elements for the transmission of a larger amount of information, and the effects of 
unwanted excitations during the process. Finally, we exploit the information-flux method to provide 
bounds to the transfer fidelity. 
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I. INTRODUCTION 

Thanks to the development of modern technologies, 
quantum information and quantum communication are 
attracting increasing interest. In this context a leading 
idea is to exploit the interaction of elements in a many- 
body register. In particular, the always-on spin-chain 
model (a one-dimensional array of two-level quantum sys¬ 
tems interacting through constant couplings) has been 
studied in detail, due to the promise of being a less chal¬ 
lenging scenario in terms of experimental realizations. It 
has been proved in theory that this system can be ex¬ 
ploited for many tasks, from the transmission of quantum 
states mm to the implementation of quantum computing 
[3] , among others. The influence that the dynamics of a 
particular element exerts on other parties during the time 
evolution of quantum many-body systems has been named 
information flux 00 . Given a certain Hamiltonian and 
the distribution of coupling strengths, by operating on 
the initial state of the system one can control the infor¬ 
mation flux and thus the dynamics. Its investigation has 
also allowed the design of protocols where the initializa¬ 
tion of the medium is not required EJE]. Unfortunately, 
many theoretical results obtained so far have not been 
demonstrated in experiments yet. This is due to the fact 
that, even if less challenging, the requirements for their 
implementations in the current setups are still demanding 
in terms of resources and control. 

With the recent progress obtained in integrated cir¬ 
cuits it has become feasible to write optical waveguides 
in glass where photons can be manipulated and interfere 
while preserving their coherence and polarization [5Ill2|. 
They have been successfully used for different applications, 
ranging from quantum walks EM to boson sampling 
fT&HTO] . These integrated waveguide circuits could also 
pave the way for experimentally demonstrating the pro¬ 
tocols put forward for spin-chain models. Even if the 
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photonic system evolves under the action of a different 
Hamiltonian, there are cases where a map to a chain of 
two-level systems is possible, and preliminary studies are 
currently in progress ESI EH. Providing an easily con¬ 
trollable setup as a proper benchmark for the theory will 
surely boost the investigation on spin-chain dynamics. In 
turn it would also be interesting, starting from the results 
already known, to generalize and modify or adapt the 
schemes to exploit the distinctive features of the bosons. 
As an indicative example, we will study in this paper 
the generalization of the state transfer to the case of a 
three-level quantum state, allowed by the larger Hilbert 
space of the single elements of our chain. 

It is important not to forget the bosonic nature of the 
systems also when we are reproducing/simulating the 
results obtained for chains of spins. This clearly has an 
effect on the efficiency of the protocols. For instance, a 
non perfect initialization of the medium used for the state 
transfer could allow more excitations to be present during 
the time evolution. As we will show in this paper, the 
behavior will thus depend on the quantum nature of the 
elements of the chain. 

After introducing in Sec. [H] the model considered in 
our investigation, we study in Sec. m the transmission 
of a general qutrit state in a bosonic chain; Sec. m 
deals with the different results obtained using a bosonic 
chain instead of a standard spin chain, when unwanted 
excitations could be present in the medium; in Sec. [V] 
we exploit the information-flux approach to provide a 
bound for the correction on transfer fidelity; finally, we 
summarize the results in Sec. EO 


II. THE SETUP: A LINEAR CHAIN OF 
COUPLED WAVEGUIDES 

In this section we would like to briefly introduce the 
experimental system we will refer to in our study. Hence¬ 
forth, we are going to discuss the transmission of infor¬ 
mation through a linear chain of N coupled sites, whose 
straightforward experimental implementation could be 
the evolution of photons in an array of N evanescently 
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coupled waveguides. The Hamiltonian that describes a first are in their ground state. Therefore, a general way 
chain of bosons with pairwise interactions reads to write the initial state reads 


N—l 

ti = E 1 + «i+l«D> (1) 

2—1 

where cq and a\ are the annihilation and creation opera¬ 
tors corresponding to channel i, and ki is the strength of 
the coupling between channels i and i + 1. We will mainly 
focus on two setups characterized by different coupling 
strengths ki, experimentally controlled by the spacing 
between waveguides: 

• the ideal case, in which couplings verify a perfect 
mirror-symmetry law: ki = J ■ \Ji ■ (N — i), where 
J is a characteristic of the energy scale, thus depend¬ 
ing on the physical implementation of the model; 

• an easier experimental configuration which cor¬ 
responds to choosing all couplings ki = J for 

[2,iV — 2] and k\ = fcjv-i = K. 

While the former is clearly interesting because it allows 
perfect state transfer (a transfer of information with unit 
fidelity) independently of the chain length [22] [23] and is 
also exploited in other protocols, the latter is a reasonable 
trade-off between a good fidelity of transmission and the 
challenges to implement the model in a realistic experi¬ 
mental setup (the values considered here for I\ are of the 
same order of J, so we are in a non-perturbative regime 
|24| ). An alternative method, exploiting local magnetic 
fields, has been proposed in |25]. The interaction time 
t, in the implementation of this model in an integrated 
waveguide circuit, is proportional to the length of the 
circuit itself. 


III. QUTRITS AND QUBITS: THE AVERAGE 
FIDELITY 

To begin with, let us analyze the approach in [T] and 
extend it to the case of qutrit transmission. As long as 
only a qubit has to be transferred and it is possible to 
initialize the chain with all the elements in their ground 
state, using a chain of bosons or fermions gives exactly the 
same results (by encoding, in the bosonic scenario, the two 
levels of the qubit as the ground and first excited state). 
The differences appear when more than one excitation 
is present in the whole chain; this could happen when 
the state to transmit requires more than two levels for its 
encoding (slightly different models have been studied in 
[26H29] ). as well as when there are unwanted excitations, 
due for instance to a non-perfect initialization of the 
medium. For the transfer of qutrits, we could encode the 
three levels as the ground and first two excited states, 
exploiting the possibility of having multiple excitations 
in the same channel. We suppose that the initial state of 
the system is separable and that all the qutrits but the 


l*(0)> = (« 10)! + P II)! + 7 |2)i) |0 ) 2 • • • 10)* , (2) 

where | 1) 1 indicates a single particle excitation in the 
first channel and |2) -! a double excitation corresponding 
to two particles in the first channel. From the definition 
of the Bloch sphere in five dimensions it follows [301 : 
a = sin 0 cos (f>e lS , /3 = sin 9 sin (f>e la and 7 = cosd, with 
the bounds 0 < 0 < f, 0 < < ■%, 0 < 5 < 2tt and 

0 < a < 27r. 

The evolved state at time t will then be 

N 

\m) =« 10)1 • • • io )*+p E 0 ' , i e_iHt i i (i ) > 1/)+ 

j ' =1 13) 

N(N+ 1)/2 V ’ 

+7 E <i"r fft i 2 ( 1 ))i/'>, 

j "=1 

where |l ( i)> = (l), | 0) 2 • ■ ■ 10 )* (|2 ( u) = |2) x | 0) 2 • • • 10 )*) 
corresponds to having one (two) photon(s) in the first 
channel and \j') (\j")) is a general state in the subspace 
of single (double) excitation. By tracing out the states 
of channels from 1 to N — 1 it is then possible to get the 
density matrix p^ ut of the qutrit in channel N at time 
t and thus the transfer fidelity. This, averaged over all 
possible initial states on a generalized Bloch sphere, reads 

F = ^2 J ^ (^in\Pout\^in) 

= 3 + 12 ^ N '^ 2 + + 1o ^ 2 ( N ), 2 ( i )\ 2 ( 4 ) 

+ ^Re[g 2(N),2(1) + fN,ig* 2 (N), 2 (l)]- 

Here, f 3 A = (j\e~’ lHt |l(i)) is the amplitude of finding a 
single photon in channel j after time t if a single excitation 
at time t = 0 was in channel 1. Similarly, 3 ^ 2 ( 1 ) = 
{i\e~ lHt \2(r^ is the amplitude of having after time t the 
double excitation corresponding to the state i (in the 
double excitation basis) if two photons were in channel 1 
at time t = 0. For the sake of completeness we observe 
that in order to compute / and g we need to apply, on the 
state the evolution is referred to, an excitation-dependent 
phase shift 

no 0 \ 

n= 0 gWigAn. 0 , (5) 

\0 0 giTr-(iV-l) J 


N being the number of channels. We remind the reader 
that, if we restrict our attention only to single particle 
excitations (qubits), the average fidelity is expressed by 


p = Re{f N , 1 ] |/iv,i| 2 1 

1 3 6 2' 


( 6 ) 


We display in Fig[l]a comparison between the fidelities 
for qubit (blue squares) and qutrit (red circles) transmis¬ 
sion through a nine-channel chain in the two scenarios 
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described in Sec. [TTJ The optimal values of t and K are 
the same for qubits and qutrits: due to the nature of 
the Hamiltonian in Eq. 0 , the transfer of two excita¬ 
tions follows the same mechanisms of the transfer of a 
single excitation, as can also be understood from the 
information-flux analysis presented later on. Constant 
experimental errors on coupling strengths normally dis¬ 
tributed and with variance a = 5% have been considered 
in the graphs. In other words, we have changed each ki 
according to ki —> ki(\ + 5), with 8 following a probability 
distribution p{8) = (l/a\/2n)e~ s / 2<T “. This corresponds 
to a scenario where the distances between waveguides are 
not exactly the ideal ones. Moreover, by introducing an 
additional Gaussian error we analyzed the effects of time 
dependence in the couplings (thus taking into account also 
the case where the waveguides are not perfectly straight). 


For a nine-channel system (with ideally equal cou¬ 
plings apart for the edge ones) we divided the evolu¬ 
tion in 100 steps each with a constant 4% Gaussian 
error plus a further 2% randomly distributed in time. 
Since fluctuations in time compensate, we obtained a 
peak in the fidelity transfer of F qu bn = (0.987 ± 0.006) 
and F qu t r it = (0.972 ± 0.012) respectively for qubits and 
qutrits, giving evidence of the system stability. It is inter¬ 
esting to notice that also the transfer of a general qutrit 
is very efficient in the regions where the fidelity for qubit 
transfer is close to 1 (i.e. for the values of the parameters 
that one should use in an experiment to obtain a good 
transmission). One could expect that this holds for higher 
dimensions as well. Hence, for the transmission of infor¬ 
mation that would require more than a single qubit for its 
encoding, it is better to exploit the larger Hilbert space 
of the chain elements and encode it in qudits, instead of 
encoding in qubits and send them one by one (that would 
require a longer total time for the transfer) or using more 
elements of the chain ED- In our analysis, we could have 
also considered the error due to photon losses. However, 
even if this is one of the main sources of error in some 
of the protocols exploiting waveguide circuits mm , in 
principle we are assuming here only straight guides. Since 
we do not have bend losses, we expect the probability to 
lose a photon to be much lower. 


To check the robustness of the protocol, we have sim¬ 
ulated the transmission with a longer chain (up to 20 
channels), finding again very similar values for the trans¬ 
fer fidelity of qubits and qutrits. It is worth mentioning 
that the second scenario (equal couplings apart from the 
edge ones) shows more robustness compared with the 
case with mirror couplings, when we consider possible 
experimental errors (we have simulated also in this case 
errors up to 5% for each coupling). This is in line with 
the results presented in [32], where it has been shown 
that more uniform couplings are more desirable. 
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FIG. 1: Comparison between the fidelities for qubit (blue 
squares) and qutrit (red circles) transmission through a 
nine-channel chain, a) The case with mirror couplings 
ki = J ■ y/i • (N — i); b) the case with equal couplings as 
a function of the edge couplings K (maximizing over 
time in the range t ■ J £ [0,20]). 


IV. QUBIT TRANSMISSION WITH EXTRA 
EXCITATIONS 

We would like now to investigate what happens to 
the transmission of qubits when the initialization of the 
medium is not perfect. In particular, we consider the 
scenario where each channel (apart from the first) has 
a probability p of not being in its ground state. As we 
want to analyze this as a case of noisy transmission, we 
consider no coherence present among the different levels of 
each channel. We restrict our study to the scenario where 
only two possible levels of these channels can be initially 
populated. In the case of a spin chain this is always true, 
but for a bosonic chain this is an approximation that 
is still good for the values of probabilities that we are 
considering, if we think about a thermal state of each 
channel. All the elements of the chain will thus start in a 
statistical mixture (1 — p) |0)(0| + p|l)(l| apart from the 
first, that will be in the 1^(0)) = a |0) 1 + |1) 1 . 
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It is now interesting to distinguish between two main 
cases. The first is the transmission of a qubit in a chain 
that potentially could transmit qutrits (or even higher 
excitations), thus including the possibility of more excita¬ 
tions in a single channel (e.g. the evolution of a photon 
qubit plus extra photons in a waveguide circuit). On the 
other hand, for the qubit transmission through proper 
spin chains, we should discard the possibility of having 
multiple excitations in the same channel. By tracing out 
the states of channels from 1 to N — 1 and averaging over 
all possible initial states on the standard Bloch sphere, 
we can evaluate the average fidelity in both the cases. It 
is a key point to stress that in theory, if we knew a-priori 
whether input errors have occurred or not, we could ap¬ 
ply different rotations on the output state: in the case 
of an extra excitation in a spin chain this corresponds 
to ili-i = |0)(0| + e I7rAr 11)(11. The transmission fidelity 
would be very close to the ideal one both for fermionic 
and bosonic systems. Unfortunately, since we can not 
know the number of excitations on the chain, we will 
always apply the rotation corresponding to qubit trans¬ 
mission TZi = 10)(0| + e ™(N-i )/2 |i((i| which in general 
deeply affects the fidelity for both systems in the case of 
multiple excitations. However, in the very particular case 
of a bosonic chain with 4 k + 1 channels ( k £ TV), or a 
fermionic chain with 4k + 3 channels (k £ N), the phase 
is revealed to be the same and thus there is no error trans¬ 
mission. As an example, in Fig|2]we compare boson and 
fermion fidelities for a nine channel chain in which each 
one of the remaining eight channels can have an extra 
input excitation with a probability p ex tra = 0.05. We can 
notice how the bosonic chain is much more convenient 
and less error sensitive for the transmission of information 
than its fermionic analogue in this case since the number 
of channels satisfies N = 4k + 1 (k £ N). 
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FIG. 2: Qubit fidelity transfer for a nine-channel chain. 
Experimental errors of ~ 5% on the couplings have been 
considered in all cases. Blue points: no extra-excitations; 
red triangles: a random extra-error excitation is averaged 
over all possible inputs in a bosonic chain; green squares: 
extra-error excitation in a fermionic chain. 



V. INFORMATION FLUX 

A quite different approach widely analyzed in jHHl 
is called information flux: it consists in working in the 
Heisenberg picture, thus focusing the attention on the 
evolution of operators, rather than the specific input state. 
Indeed, we would like to understand how this information 
flux is related to the fidelity and in particular if it could 
be used to estimate the efficiency of the transmission in 
more general cases than those studied so far. 

In the Heisenberg picture the evolution of an op¬ 
erator is regulated by Oflt) = W(t)OiU(t), where 
U(t) = Texp[—(i/H) f i-L{t')dt']. In the case of time in¬ 
dependent Hamiltonian, the expression is simplified in 
Oi{t) = Oie~Ti nt , which can be easily expanded as 

6i(t) = o l + A] + [HA]] + ■■■ (?) 

As it has already been shown BH3, recalling (JT]) and 
computing a series of commutators, it is possible to write 
the evolved operator Oflt) in terms of a set of time- 
dependent parameters C(t) and a set of operators Qj (all 
considered at time to = 0), arising from the commutators. 

Hence, we can now evaluate the time evolution of cre¬ 
ation and annihilation operators in the case of bosons 
evolving with the Hamiltonian Q. What we find out is 
that the set of commutators is composed respectively by 
all the creation and annihilation operators themselves. As 
a result, also the evolution can be expressed in the same 
terms: 


N N 

a\(t) =^2cn^ j) {t)d](0 ) , di{t) = ^2AN {i<j) {t)dj{Q). 
j =i 3=1 

(8) 

Given the evolution of a certain operator acting on 
the last channel of the chain, we define information 
flux as the coefficient that sets the weight of the same 
operator applied at time zero on the first channel: 

X°^ a = AJ\T(ClZ)( Nt iy We highlight that in this case 
C7£(t)(jj) = A\f(t)yj) V i,j,t for obvious symmetry 
reasons: for the sake of simplicity we will denote from 
now on as Cyy both CIZ^ ^ and AAfyjy 

The evolution of the observable related to the number 
of photons is 

N 

d\(t)di(t) = y ClZ(ij)(t)AN(i,k){t)a\a k (9) 

3,k=l 

and the corresponding information flux 

xy = cn (Ntl) AAf (Nil) . ( 10 ) 

In Fig j3] we display the behavior of the information 
flux T'j^" in a few relevant cases that we have already 
examined: the one with mirror couplings and with equal 
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FIG. 3: Recalling Figjljwe plot the information flux Xy 0 
for a nine-channel system in the case of mirror couplings 
(a) and for equal couplings (b) (maximizing over time in 
the range t ■ J £ [0, 20]). Experimental errors of 5% on 
the couplings have been considered. 


couplings (experimental errors of the order of 5% on the 
couplings are taken into account). 

We remark that if we encode a qubit (qutrit) through 
the 10), |1) (| 2)) basis the number of photons in a channel 
is strictly related with the fidelity: we show in Fig|4] a 
comparison between information flux and fidelity for qubit 
and qutrits in a chain with equal couplings within a 5% 
error. In addition, in the next section we would like to 
better understand how the ajya.y (f) operator that acts on 
the last channel at time t is related with the information 
flux by considering various possible initial input states 
for channels 2 • • • N. In particular, we will focus on Fock 
states, coherent states and a more general undefined state, 
though always considering the states in every channel pure 
and separable with respect to each other. 
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FIG. 4: Information flux 2jy“ (green squares) together 
with the transmission fidelity for qubits (blue circles) and 
qutrits (red triangles) in a linear chain with 9 channels 
(maximizing over time in the range t ■ J £ [0, 20]). 


A. Fock states 

Let us call |SRztv( 0)) the initial state of the system, 
being indeed a separable state where we have encoded 
our qubit in the first channel while every channel * > 2 is 
initially in the Fock state |z(0)). We would then have 

n n(t) = IN{ty\a ! N a N (t)\'f> i N ( 0 )^ 

= Y. CN,jCN,k (^IN (0)1 a]0/^1^777(0)^ 

j,k 

= E c 'iW°) (11) 

3 

N 

= li a nM + Y.C 2 N,M 0 ). 

7=2 

The number of output photons in the last channel is thus 
given by the initial number of photons in the first channel, 
weighted by the information flux X y a , with a correction 
depending on the number of photons initially present in 
the other channels. 


B. Coherent states 

This time we suppose the system initially in a separable 
state where |ay) = c n)j is the initial 
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coherent state in channel j > 2. We would then have 

^ CN,jCN,k ^U'/iv(0)|ata fe |^'/j V (0)^ = E CN,jCN,kOt*jak 

j,k j,k 

N 

= C% tl ni(0) + yy ^N,j\ a j\ 2 + 2 yy i?e[a*afc]Ujv,jCjv,fc+ 

j=2 j>k>l 

+ 2 yy iZe[a*]\/ni(0)C' J v,jC'jv,i- 

j¥1 

( 12 ) 

By manipulating these terms (as shown in the Appendix), 
we can find an upper bound for their weight, thus con¬ 
cluding that ^/Ar( 0 )|a] v ajv(t)|^ , 7 Ar( 0 )^ = ni(0)C^ 1 +r. 
Here T represents the corrections which are bounded by 

r<\a max \ 2 (i-c 2 NA )(N-i)+ 

+ 2^nA0)C NA \a max \VW^l^(l - C 2 N1 ), (13) 

where we point out that C^ 1 = °. 


one of the two parentheses is the result of the product 
i (7(0)|a|i(0)) i on the initial state of channel i. Since the 
state |i(0)) is normalized in the Fock basis, this quantity 
is upper-bounded by 1, and it is actually equal to one only 
in the case of a coherent state, which therefore reveals to 
be the worst case. This means that the second term can 
be bounded by: 

(^7iv(0)|a]a fc |^7jv(0)) | j^k < E CN,jCw,ky/njVnk 

a 8 ) 

Following exactly the same procedure of the previous 
paragraph, we can finally conclude that in the general 

case (^iN{0)\a l N a N (t)\'i’ IN {0)^ = C 2 Nl -ni(0) + T. Here, 
r represents the corrections which are bounded by 

T < rij(0) mox (l - C% a )(N - 1)+ 

2^nA0)C N ,i^/n J (0) max VN^l^/(l - C 2 N1 ), 

where again C^ 1 = X^ a . 


C. More general states 


Let us call |H> 7 tv( 0 )) the initial state (at time t = 0): 

N 

IH>7tv( 0)) = n^(O) 10)., + A(0) |1), + 7i(0) |2>. + •••), 

i =1 

(14) 

where i represents the channel and a(t), /3(t), 7 (f), • • • 
are the coefficients at time t respectively of the vacuum 
component, one particle component, etc. 

Now, recalling the general expression 

^7iv ( 0 )| 0 ^ 077(^1 ^ 7 ^( 0 ) = 

= yy cnjCn^ (*i N m 4 ak\*iN(o)), 

j,k 

it is possible to split it into two parts 

{^(0)|ota j |$7 JV (0 ))=n j ( 0) (16) 

and 

^4 , 7Ar(0)|a]afc|'F7Ar(0)^) | j^k = 

= (M 0)4(0) + \/2 7fc (0)/?£(0) + 734(0)7^(0) + •••)• 

• (aj(0)(3* (0) + yfifr (0)7* (0) + T3 7i (0)y*(0) + •••)• 

(17) 

While the first term is easy to calculate, being actually 
very similar to the case of coherent states, instead, in 
order to maximize the second term, we observe that each 


VI. CONCLUSIONS 

The rapid and striking progress obtained in the field 
of integrated circuits has allowed the experimental im¬ 
plementation of several quantum information protocols. 
Another promising direction for exploiting these setups 
would be the adaptation of schemes that have been already 
proposed in literature for spin chains. Considering their 
wide range of applications, from the faithful transmission 
of information to its manipulation and control, including 
the possibility of generating quantum entanglement, this 
will clearly pave the way to further progress. However, 
we have to keep in mind that, even if the map from a 
fermionic to a bosonic chain is exact in some particular 
scenario, this is not always the case. One can also find 
advantages in using these adapted protocols with respect 
to the standard spin-chain schemes. For instance, here 
we have investigated the possibility of transmitting three- 
level quantum states by encoding them in the multiple 
channel-excitations of a bosonic chain, with a fidelity 
very close to the case of qubits in spin chains. Moreover, 
we have analyzed the effect on the transfer fidelity of 
unwanted excitations, that could happen in a realistic 
experimental implementation. We have highlighted how 
the cases of bosonic and fermionic chains give very similar 
results, but they could strongly differ if the number of 
channels satisfies particular conditions. Finally, an analy¬ 
sis from a different viewpoint, namely an approach based 
on the information flux, has allowed us to obtain bounds 
to the fidelity in more general scenarios. 
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Appendix A: Coherent states 


Hereafter we report the calculations of the correction term for information flux in the case of coherent states in 
input. 


N 

E 

3 =2 


c 2 3 4 5 6 

U N,j 


(*j | < \Oj 


(1 — Cjv ,i) 


(Al) 




2 ^2 R e [ a j a k\CN,jCN,k < |e^max| 2 E 

r --— . - a/(1 — Cat,i) 

< Kax| 2 • ^(1 - U^ 1 )V / F^T . V . (V - 2) 

= |a ma:c | 2 (l-C'^ 1 )(^-2) 


(A2) 


where, since y~V | C| 2 = 1, the maximum of the sum i |Cjv,j| is obtained in the case in which the coefficients 

are all the same and equal to |Cjv,j| = y — n-^ 1 • addition, the second inequality holds since the second sum 
covers N — 2 terms, while the first sum is over TV — 1 addends. 


2 \J ni(0) Re[a*]C N jC N 'i = 2y / m(0)(7jv,i ■ '^2 a *i 

i 

G 2 \pfV\ (0)Cjv,i |o; max | El^dl (A3) 

j¥= i 

< 2 V / ^(0yC J v, 1 |a ma;c |v / TV^iy(l - 

where again, since CV.j | 2 = 1, the maximum of the sum I Civ, j I is obtained in the case in which the 

coefficients are all the same and equal to Cfy.j = 
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